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Abstract 

Data from experimental observations of certain brain functions, such as impulse responses of cortex to electric shocks (av¬ 
erage evoked potentials), related with the coexistence of pulses and wave modes in the dynamics of a class of neurological 
processes (Freeman K-sets), present functional distribution reproducing the Bessel function behavior. This suggests the possi¬ 
bility to replace ordinary differential equations, typically used to model data bases concerning such processes, with couples of 
damped/amplified oscillators which provide time dependent representation of spherical Bessel equation. The root loci of poles 
and zeros of the equation solutions are shown to conform to solutions of K-sets. One advantage of the present formalism is 
that some light is shed on the problem of filling the gap between the behavior at cellular level and the macroscopic dynamics 
involved in the traffic between the brain and its environment. Breakdown of time-reversal symmetry in each of the (damped 
and amplified) oscillators is related with the cortex thermodynamic features. This is proposed to be a possible mechanism to 
deduce lifetime of recorded memory. 
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I. Introduction 

There is an essential problem in the study of brain function that even today, after so many years since Karl Lashley posed 
his dilemma, still waits for a solution. As recalled many times in the literature, in the mid 1940s he wrote [J|: 

” .... Here is the dilemma. Nerve impulses are transmitted ...from cell to cell through definite intercellular connections. Yet, 
all behavior seems to be determined by masses of excitation...within general fields of activity, without regard to particular nerve 
cells... What sort of nervous organization might be capable of responding to a pattern of excitation without limited specialized 
path of conduction? The problem is almost universal in the activity of the nervous system.” 

The successes of neuroscience in the study of the structural and biochemical properties of neurons, glia cells, and all the 
biological units and cellular structures in the brain have not yet filled the gap between the behavior quite well understood at 
cellular level and the macroscopic dynamics involved in the traffic between the brain and the world around it. Much work is 
still to be done in order to have better and better insight in a matter where the connection between phenomenological data and 
a reliable theoretical representation is still facing extreme difhculties. Within the limits discussed below, this is the problem on 
which we focus our attention in this paper. 

In order to analyze it, we consider data from experimental observations made by fitting measurements of brain function 
with solutions to ordinary differential equations (ODE) (Freeman K-sets) @|. From the observations it emerges that data 
from neurological processes, such as impulse responses of cortex to electric shocks (average evoked potentials (AEP)), present 
functional distribution in plots of amplitude (in fiVolts) vs time (in msec) as shown in (Fig. 1) Q] reproducible by the Bessel 
function behavior. Since Bessel equations have a time dependent representation in terms of couples of equations, each couple 
containing one equation for a damped oscillator, the other one for an amplified oscillator, both with time dependent frequency, 
we are motivated to study the observed neural processes in terms of these oscillators, which is also consistent with the fact 


2 


that consideration of dynamics in the continuum of the wave mode cannot be omitted, as it appears from the experimental 
observations we start with, as well as from other data and analysis Ml (see also Appendix A) showing that pulses and wave 
modes coexist simultaneously in the dynamics of neural populations, whose understanding thus requires to consider both of 
them. 

The behavior of each parametric oscillator emerges from an interactive population of excitatory cells coupled with an 
interactive population of inhibitory cells forming a negative feedback loop, in which the two populations are in quadrature 
having the same frequency and decay rate, with on average the excitatory oscillation leading the inhibitory oscillation by 
7r/2 radians (cf. Appendix A.3)^. The damped oscillator has low feedback gain; the amplified oscillator has high feedback 
gain. We obtain, as a result of our study, the agreement between the dissipative dynamics expressed by each oscillator in 
the couples and the conclusions on the root loci and the evolution towards limit cycles presented in @] (see Section IV and 
Appendix B). Our analysis in terms of the parametric oscillator couples is also consistent with temporal frequency modulation 
as it is derived from observations and widely used in the literature . In addition to this, we also obtain a link 

with the dissipative many-body model of brain Si, since the coimles of damped/amplified oscillators is known to have a 
representation in terms of correlation modes in many-body systems [8-[ll||. Each couple forms an operator and its double, one 
running forward in time in the brain, the other backward in time in the mind. 

Our analysis is limited to the data mentioned above which are well fitted by the Bessel function. Perhaps, a plausible 
macroscopic physical justihcation of the observed similarity of evoked potentials with Bessel functions is likely due to the 
2-dimesionality of the distribution of the excitatory and inhibitory neurons. On the other hand, the occurrence of Bessel 
functions might be not so surprising since in physics Bessel functions play a crucial role in theoretical descriptions of (resonant) 
pulses in cross sections of particle scattering as well as of wave propagation processes, which also suggests that the observed 
brain functions have, indeed, a lot to do, simultaneously, with pulse modes and wave modes, as already mentioned. The 
representation in terms of couples of damped/amplified oscillators also allows the interplay of linearity/nonlinearity dynamics 
features [l^. Such an interplay is of crucial importance in brain studies (see our comments in closing Section II). Linearity 
allows us to consider the summations in the fitting functions, as we will see in the following Sections, and it has been one of the 
motivations to adopt the analysis in terms of ODE in Q]. On the other hand, the biophysical nonlinearity at the microscopic 
level appears at the macroscopic level as the root loci^ revealing the capability by neural activity to provide the transition 
energy necessary to sustain the phase transition. Such a nonlinearity is present in our modeling since the damped/amplified 
oscillators are actually nonlinearly coupled as it appears considering their dynamical properties fia . [a. We will comment more 
on this point in Section II. 

In this paper, we do not consider plots of evoked potential data showing functional distributions different from Bessel 
functions. However, since any function can be expanded in a series of Bessel functions (see e.g. [IM), the representation 
in terms of damped/amplified oscillators and our results can be extended also to such cases, within quite general boundary 
conditions and provided a convenient analysis is carried on, which we plan to do in a future work. 

The functional activity of the brain and the laboratory observations at the cellular (neuronal) level represent the manifestation 
of the dynamics underlying at the elementary constituent level. The degrees of freedom of such a dynamics are associated with 
the elementary components and modes of the molecules present in the system, such as the electrical dipole vibrational modes 
of molecules of biological interest and of the water molecules constituting the bath in which they are embedded. These are 
of course quantum degrees of freedom. Their dynamics is studied in the dissipative quantum model of brain m , which deals 
indeed with the quantization of the damped/amplified couples of oscillators. We do not consider the dissipative model in this 
paper. In any case, we stress that neurons and glia cells are considered to be classical objects in this paper and in the dissipative 
model. 

The plan of the paper is the following. In Section II we outline the derivation of the time dependent representation of the 
Bessel equation in terms of the damped/amplified couples of oscillators. We present in the Appendices A the details of the 
empirical setting in which the data analyzed in this paper were derived. Section III is devoted to the analysis of the data 
in terms of the couple of damped/amplified oscillator equations. Applications and explicit examples of computation in our 
formalism are discussed in Section IV, with formal details given in the Appendix B. In Section V, a hierarchy of lifetimes 
for k-mode (k is a wave number) components is considered. Section VI is devoted to conclusions. For the sake of brevity, in 
this paper, we do not discuss features such as scale free behaviors in brain functioning, which may be analyzed in in terms of 
coherent state dynamics in the dissipative model. For a discussion on this last point see Refs. [itI. H^. 


II. Bessel equations and the damped/amplified oscillator couples 

A conclusion that can be drawn from the laboratory data and their analysis, as summarized in the Appendix A where the 
empirical setting and the scenario in which the data were collected are presented, is that the (Wnamics of neural populations 
cannot be understood from either waves or pulses alone. As well known in the literature 0, y, [3, pulses and wave modes 
indeed coexist simultaneously in circular causality. Each determines the other. Laboratory observations are made by batch 
processing, with categorization after each shot. What has been found Q] is that, in the successive brief time windows accessible 


1 


The phase lead may be varied from average by the two kinds of positive feedback in the populations if they are not equal to the negative 
feedback gain. 
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FIG. 1: AEP and fitted curve for impulse response of the prepyriform cortex in the closed loop state. Picture adapted from 
Ref. [2], Section II, Page 111. 


to observations, the dynamics of dendritic integration and axonal transmission by neural populations combine the two modes. 
On the one hand, sometimes the paradigm prevails which uses only data from the pulse mode. It omits consideration of 
dynamics in the continuum of the wave mode and focuses on the point processes revealed by action potentials. On the other 
hand, one cannot omits consideration of dynamics in the continuum of the wave mode. The above-mentioned laboratory 
observation that functions representing neurobiological processes can be fitted by Bessel functions is significant for opening to 
investigation the properties of neural collectives. We consider thus the spherical Bessel equation of order n [l3-[l^: 


-f 22—-b [z^ - n{n+l)]Mr, 


( 1 ) 


with n being zero or an integer. Equation 0 admits as particular solutions Bessel functions of the first and second kind. 
Superposition of both is called a Hankel function [iMl. 

Note that a transformation of the type n —(n+1) leaves Equation ([T]) invariant so that both and are solutions 

of the same equation. This degeneracy can be lifted by introducing the parametrization as done in [2^: qn,i ~ Mn ■ 

_ t 

z = z{t) = tnC , where and a„ are arbitrary parameters and t denotes the time variable. e„ is dimesionless and 
has the dimension of time. By setting n(n + 1) = 1(1 -\- 1), one obtains two differential equations for I = —(n -|- 1) and I = n, 
respectively, where the symmetry under the transformation n —>■ —(n -|- 1) is broken [^. These differential equations can be 
reduced to the parametric damped/amplified oscillator equations: 


^+U + ujl(t)^ = 0 , (2) 

C-Li + Lol(t)C = 0 , (3) 


where the dots represent derivatives with respect to time, L is a parameter intrinsic to the system (e.g. a biological, physiological 
or behavioral parameter) and the shorthand notation is: ^ = q„^_(ri+i) , C = Qn,n , T = ^2'^^ • time-dependent frequencies 

Lt _ Z- 

arising after introducing the above parametrization is uJn{t) = ^n+i = uj^e . Note that loq = = kc = a;o,fc, with k 

the wave-number k — 2 -k jX. The k dependence of L can be expressed as Ls, = This is done in such a way that e„, Ls, 

and cuo,fe, which have dimension of , are independent of n (fixed time-invariant frequencies). Note that the transformation 
n —(n -b 1) leads to solutions (corresponding to M_(„_(.i)) which can be obtained by time-reversal t —^ —t in u!n(t). We will 
omit the k index for simplicity whenever no misunderstanding arises. 

Eqs . IZl l and © are recognized to be nothing else than the couple of equations at the basis of the dissipative model of 
brain [1^ with time-dependent frequency iij(t). They can be reduced to the single parametric oscillator 

r„-b fl^(t)r„ = 0 , (4) 


1 _ Lt L t 

where = -^r^e 2 , 2 and tin is the common frequency of oscillation: 

tln(t)=[ujl(t)-^'\^ >Q , 


(5) 
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so that the functions qn^n and are “harmonically conjugated”. Proof of the mechanism is that the signals are observed 

to be in quadrature; the frequency and decay rate are always identical, and the phase of the inhibitory signal lags a quarter 
cycle behind the excitatory signal. Note that i}{t) becomes time-independent as t —^ oo and it is fully specified by characteristic 
parameters of the system^. Q is assumed to be real in order to avoid the overdamping regime. 

In conclusion, the Bessel equation ((T|, for each n, has a time dependent representation in terms of the Eqs. ^ and (0 of 
the couple of damped/amplified oscillators, provided the above parametrization is adopted. 

We observe that in the dissipative quantum model of brain Eqs. and m describe electric dipole vibrational modes, denoted 
by flfc and at, respectively, called dipole wave quanta (DWQ), with the suffix k denoting the wave-number introduced above. 
These are the Nambu-Godstone (NG) boson quanta dynamically generated as a result of the breakdown of the rotational 
symmetry of the molecule electrical dipoles & ij, [2j| . Their condensation in the ground state of the system leads to the 
formation of ordered (coherent) patterns [23-[^|. The brain couples to its environment (referred to as its Double 
through the perception of external stimuli. These produce the symmetry breakdown and trigger the condensation process. This 
is constrained by the requirement J\fa^, — Afaj, = 0, with A/’aj, and A/a^ denoting the number of au and Ofc modes, respectively, so 
that the ground state has the lowest energy. Such a constrain does not fix, however, the number of fc-modes, for each k, in the 
ground state and we can have infinitely many lowest energy states with distinct combinations of different momentum modes, 
characterized by the “order parameter” M = {A/'aj, = A/a^,, V/c}. The A/’-set thus labels the system ground states (referred to 
as the the memory states in the dissipative model). 

We close this Section with two comments. The first one is about the breakdown of time-reversal symmetry introduced by 
Eqs. © and ([S]). This is intrinsic to the dissipative/amplified character of Eq. ([2| and Eq. ([21). Breakdown of time-reversal 
symmetry introduces the “arrow of time” (in each of the damped and amplified oscillators), which, as discussed in Refs. [23-1^. 
can be related to the “thermodynamical” and the “psychological” arrows of time. Among the many existing papers dealing 
with nonequilibrium thermodynamics and time-reversal symmetry breakdown in neuronal systems we only quote here few of 
them, see [25l - [^ . We do not insist further on this topic in this paper. 

The other comment concerns the fact that the system of equations and © allows to us to consider the summations 

in the fitting functions, as we will see in the following Section, cf. Eqs. m and ©• Linearity, which, as recalled in the 
Introduction, has been one of the motivations to adopt the analysis in terms of ODE in Q| , is thus implicit also in the present 
treatment. For observational analysis is highly convenient that the generic equation for the impulse response in time and 
space is the sum of the responses in the two modes of observation (pulse and wave mode). This motivates the choice of the 
fitting functions used in Ref. and considered in the next Section, see Eq. (|21) (and (|S])). On the other hand, biophysical 
nonlinearity appearing at the macroscopic level as the root loci @| reveals the capability by neural activity to provide energy 
necessary to sustain phase transition resembling a Hopf bifurcation. It is well known that nonlinearity is crucial in brain, 
especially in the pulse-wave transit asymmetric sigmoid function [^, [^. Such a nonlinearity is also present in our modeling 
since the damped/amplified oscillators are known to be nonlinearly coupled [12], L being the dissipation/loading parameter in 
the combined damped/amplified system dynamics, as it appears also from the entangled two modes quantum coherent state 
they generate at quantum level (the interested reader is referred to [12, [2l|| for details). Moreover, the root loci studied in 
Ref. [3|, and in the literature therein quoted, also appear in the present treatment, as we will show in the following Sections. 
It is also interesting that the time independent value wo of the frequency is reached for higher values of n, namely in the 

limit of the maximal coupling between the damped and the amplified oscillators (the whole system (^, () is a closed system 
in that limit of large n. Considering that a brain is a thermodynamic system with a fuzzy number of ports, n as an integer 
must be very large indeed.). This is the limit of stationarity: observations show that most cortical impulse responses (evoked 
potentials from single shock electrical stimulation) lasted only 1.5 to 2.5 cycles, and those few lasting longer had fixed stationary 
frequency. The linear range is crucial for the test of superposition (which allows the summation in the fitting functions © and 
(|5|)) since there the pulse density can be represented as a function of wave density through the G(p) and/or G(v) functions (cf. 
Eqs. (IA3I) and (IA2I) and Ref. [2]). In order to use time ensemble averaging to estimate pulse density in the linear range, one 
also assumes ergodicity, which can be done also in the present case. Thus, linearity, stationarity and ergodicity, combined with 
nonlinearity, are compatible with the damped/amplified parametric oscillator treatment. 

Our next task is to show how the formalism concretely works in explicit examples. 


III. Fitting functions and the Bessel equation 

In this section we show that neuronal processes described in Ref. [2 can be studied by using the damped/amplified oscillator 
equations along the lines discussed above leading to results in agreement with those obtained in Qj. In Section IV and V 
we then present explicit examples of computations (see also the Appendix B) and discuss the lifetime of correlation modes 
(fc-modes), respectively, in order to show how the present scheme concretely works. 

The fitting functions for the laboratory observations used in [2| (see in particular Section 2.5.3) to describe the oscillatory 


^ Gamma bursts of oscillation in the ECoG lasting long enough (> 3 to 5 cycles) commonly show temporal frequency modulation (smooth 
decrease of 10% of the mean, but rarely an increase). The impulse responses (average evoked potentials) has fixed time-invariant 
frequencies. This is equivalent to time t going to infinity. 
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processes resulting from connections between sets of neurons are generally of the form 


v{t) = Y^Aje * , (6) 

3 


where Aj, pj € C with positive coefficients (i.e. Re{pj) > 0 and Im{pj) > 0). When pj is real, Aj is also real. Otherwise, 
Pj and Aj always come in conjugate pairs to yield damping sin or cos terms of the type sin * with ipj 

being a constant phase term in Aj. For example, using the Euler identity 2sin0 = — ie*®, with ‘angle’ 9, the term 

sin(rj t + (pj) * can be rewritten as Aje~^^ * + * where Aj = '2^/2, pj = 7 j and similarly for a damping 

cos term. We observe that Aj and pj actually depend on the wave-number kj. In the field theory frame adopted by us, kj, 
for each j, varies over a continuum set of values. We will use the simplified notation Aj = At- and pj = pk^ whenever no 
misunderstanding arises. The Laplace transform of Eq. m is 


V{s) 


n"(^+Pi) ’ 


(7) 


which has the zeros at the at Zi and all its poles lying in the left-hand side of the complex s-plane. 

Each term in the summation in Eq. corresponds to a DWQ ak^ mode of the dissipative quantum model of brain. This 
is suggested by the fact that, as we will explicitly show, each term in Eq. m satisfies Eq. ©■ Hence, the modeled output 
potential is described by a linear superposition of at - modes which represent solutions of the n-order Bessel equation. In this 
way, a particular stimulus will excite an A/'-set which is made up of different A/"aj, components f^ . 

Moreover, one should remark that, for each j, the complete basis of functions, of which * ig one element, is {e~^^ *, 

This means that one cannot avoid to consider also the element other words, mathematical completeness requires that 

we need to “double” the degrees of freedom of the system by introducing the other functions of the basis, and thus we 

also need to consider 


v{t) ’ ( 8 ) 

j 

where again Aj, pj are complex numbers with positive coefficients. These terms are solutions of Eq. Q associated to the 
ak modes. This shows the consistency between the mathematical necessity of considering the complete functional basis 
and the derivation in Section II of the couple of damped and amplified oscillator equations. It is worth noting 
that the Laplace transform of Eq. m is 


yi;{^-p3) ’ 

(9) 

which has poles lying in the right hand side of the complex s-space. Restoring the k subscript notation, 
Aj / 0, in Eq. ([2| gives: 

use of Aj(t) = Aj 

plj - Lk^Pkj + ujl^kj (f) = 0 , 

(10) 

so that 


Lk 

Pkj = 2 ^ ,n{^) 5 

(11) 

for real coefficients. For complex coefficients we obtain: 


7fc,- - rlj + 2*(ifej Tfej) - Lkj (jkj + *r kj )+ujl_kj{t) = 0 . 

(12) 

After separating real and imaginary parts, 


Lkj = 2'ykj , 

(13) 

= \pkj\^ , 

(14) 

and rearranging: 


= -Y ’ - 7fe,] "" = . 

(15) 

So that Eq. (USD agrees with Eq. (Hill and in both, ^kj,n{t) is given by Eq. ([5|l. 

By doubling the degrees of freedom, as required by the completeness of the basis {e~^^ and introducing terms from 

the summation in Eq. ([8| as solutions of Eq. Q, one obtains the same equations for the coefficients as in Eqs. Hill and H5II. 
with the same common frequency 0.kj,n{t), which thus establishes the (resonating) link between the damped/amplified systems 
for each fe-mode and a particular n (cf. Section II). We are now ready to work out explicit examples. 
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IV. Examples of computations and limit cycle 

The motivation of the present Section is to show that the representation in terms of parametric damped/amplihed oscillator 
couples leads to results which are in agreement with the ones obtained by use of ODE [^. To this aim, we study two explicit 
examples from Ref. Q], Section 2.5.3. We then discuss briefly the existence of limit cycle. Since computations are straightforward 
we summarize them in the Appendix B. 

The first example we consider is the one of the average evoked potential from the prepyriform cortex of a cat. It is fitted by 
means of a non-linear regression to yield the following equation (KII set) [§]: 

v{t) = 14i[sin(rfcjt-|-(pfcJe“^'“i* -sin((pfcje“'’''i*]-|- 

-I-142[sin(rfe2t-I-- sin((^fe2)e“'’'“2‘] , (16) 

where each term corresponds to a particular value of momentum kj, j — 1, 2, associated to each mode. 14^ 2 are constants, the 
T’s and 7 ’s correspond to the frequency and decay rate of each component respectively, /I’s are the rates of change of the initial 
peak and their values can be found experimentally. The interesting aspect of this particular example is that one can postulate 
that the /3 value is in fact slightly different for the dominant and the subsidiary component (1 and 2 , respectively) and each of 
these components is a solution for the same fc-value, as it is the damping sinusoidal term. Details of the computation in the 
present formalism are reported in the Appendix B.l. 

One can perform the analysis of Eq. (HU by plotting its zeros and poles in the complex plane as done in Fig. 2. We observe 
that a slight change in the values produces changes in the shape of the fitting function (See Fig. 3). This, together with the 
role of P in matching the growing and decaying solutions (cf. Appendix B.l), means that the brain might use the P parameter 
to discriminate between two different, though similar, behaviors or perceptions. 

In the Appendix B.l we also study the addition of two different fc-modes with an output of the form of Eqs. (IB12I) . We find 
that if we have a linear superposition of two fc-modes of the form of Eqs. HU, where each satishes Eq. © with a different 
value of P, the addition of the two modes will present 4 complex poles and 2 complex zeros in s-space. In general, each extra 
mode of the same type added to the output function increases the number of complex poles and complex zeros by two when 
the output is represented in s-space. 

As a second example, we consider the modeling of the synapses of axons in the primary olfactory nerve (PON) (cf. also 
[^, Section 5.2.2). These synapses are of particular interest. The model that represents them has to account for the time delay 
between input and output. 

We compute the root locus plot and can reinterpret the transfer functions in t-space as a linear superposition of the nonlinear 
fc-modes. These coefficients will be given by the appropriate parameters depending on whether they are real or complex. 
Moreover, they can be determined by the position of the poles of the transfer function in s-space which are gain-dependent 
(Fig. 4). 

Another interesting feature is that the doubling of the degrees of freedom produces a mirror image of the poles in s-space 
with respect to the imaginary axis. The pole distribution and its variations with gain is critical when determining the chaotic 
dynamics of the neurological responses [^. An important point has to do with the “crossing ”of the poles from the left to the 
right hand side in s-space @|. 

The connections between the left and right hand side of the locus plots (i.e. damped and amplified solutions) seem to be 
deeply related with the chaotic dynamics that is observed in the brain [^. [ 3 ^. 

In the Appendix B.2 we reproduce one of the examples in Ql and link it to the present formalism, making thus explicit its 
connection with the analysis of [ 2 ]. 

The root locus plot of the transfer function D{s) given by Eq. (IB19I) is represented in Fig. 5 and one can see how the larger 
poles cross the imaginary axis from the left to the right as decreases, which is a sign of an unstable limit cycle. The limit 
cycle will be stable at frequencies where the poles are purely imaginary (« 160 rad/sec). Within our framework, the crossing 
of the imaginary axis means the transition from a decaying solution to an amplifying one as the gain increases. 

The existence of a stable limit cycle means that the rate of oscillation of the pulses (Ffc^. ^„) changes with time at a fixed 
frequency. This agrees with our model and identifies the limit cycle frequency with Q.(t) which in fact, links the damped and 
amplified solutions. Therefore, the stable equilibrium state (constant pulse frequency) corresponds to fully coupled solutions 
(n —^ 00 ) since fl(t) becomes time independent in this limit (see Section II and next Section). In addition, the stable limit 
cycle occurs for solutions with = 0 . 


V. Correlation modes and their lifetime 


Our discussion has been focused on damping and amplifying terms in the fitting functions and on the root loci of zeros and 
poles. Now we want to discuss about lifetime and spatial range of correlation modes. 

We start by recalling that the common frequency n(t), Eq. ([U, has been assumed to be real at any time t in order to avoid 
overdamping. Inspection of Eq. shows that such a reality condition of D(t) imposes an upper bound, say Tk^,n, in the time 
within which a process can take place given by 



( 17 ) 
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FIG. 2: Poles (x) and zeros (•) in s-space from Eq. (lB6ll . 



FIG. 3: Fitting functions with same (thick/continuous) or different (dashed) /3 values. 


L ^^3 

Eq. HZl) imposes a bound on the admissible values of k\ using oJo^kj = kc, we have k > = k(n, t). This can be 

equivalently expressed in terms of a cut-off on the linear range A oc 1/fc over which the fc-mode can propagate (wave-lengths 
greater than A oc l/fc(n, t) for any n and t are excluded). It it is useful to introduce the quantity \kj,n{t) [^. 1^: 

In^sinhj^Tfe.,™) -ln('sinh2^(Tfe^.„-t)W4^t 

r^kj,n\t) — ^ , (18) 

so that Afc^.,n(0) = 0 for any k and Kk-,n{t) oo for t Tkj,n for any n. Then, can be written as; 

fIfc.,„(Afc,.n(t)) = . (19) 

Therefore, i3kj,ni-^kj,niTkj,n)) = 0 so that Kk-^n oc Tk^^n = lifetime of each fc-mode. In this way we see that each fc-mode 
“lives” with a proper lifetime rk-,n, so that the mode is born when Tk-,n is zero and it dies for Tk^^n —^ oo [2^. This is illustrated 
in Fig. 6 for growing values of k and fixed n = 1. The figure, where different Ak^,n are plotted vs t, shows that the lifetime 













FIG. 4: Poles (x) in root loci for 0 < JFe < 1 from Eq. (IB15II and zeros (•) from Eq. (IB17II and its “Double”. The poles at 
s = ±2300 and the two zeros at s = ±2300 are not shown. 


jtu 



FIG. 5: Root locus plot of Equation (IB19I) . 


increases as k increases. According to this model, the modes that do not satisfy the fl(t) reality condition are decayed modes 
i^kj,n = 0), which agrees with the requirement in @] that the damping coefficients describing the voltage output must be real 
and positive (cf. Eq. dill) '). Since smaller k implies larger A, shorter lifetime modes condense over larger domain. Vice-versa, 
larger k will condense on smaller domains, which, however, live longer. 

The sinusoidal terms in the fitting functions are modulated by the time-dependent frequency Q{t) so that, when the particular 
mode decays (fl(t) = 0), the term becomes purely damping/amplifying. On the other hand, the damping/amplifying coefficients 
are always expressed in terms of L, which we have seen is determined by intrinsic parameters of the system. In addition, the 
frequency D(t) becomes time independent when n —>■ oo. In such a limit, the sinusoidal terms of the fitting functions describe 
the system dynamical regime where energy is conserved and the modes auj fully couple to the modes a/c^.. As already observed, 
n represents the number of links between a and a and we see that, as the number of links increases, the lifetimes of the fc-modes 
are longer (Fig. 7). 

Finally, the cooperative role of n and k in relation to the lifetimes is represented in Figs. 8 and 9 for “low” and “high” values 
of k In relation to fco, respectively, for growing values of n (the units were chosen so that c = 1). 
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FIG. 6: ^kj,i vs t for different values of k and fixed n = 1. 



FIG. 7: A 3 ,n vs t for different values of n and fixed fe = 3 in units of c = 1. 


VI. Conclusion 

The lesson we 
only progress in 
them both, since 

hand, it is the nonlinearity in the pulse-to-wave conversion that explains the properties of stability and instability relating to 
criticality and phase transitions and the variable couplings among neurons that explain reinforcement learning in populations. 
However, most of the time most neurons in most populations are functioning in domains of superposition, whether governed 
by a point attractor in a basal chaotic phase or by a limit cycle attractor in an orderly active phase. An observer faced with 
understanding such a complex dynamical system is well advised to seek for a domain in which the dynamics is linear, stationary 
and Gaussian, as commented upon in previous Sections (and in Appendix A). The existence of domains of superposition in 
the cortical dynamic state space is indeed exceedingly important. Analytic solutions of the equations for multiple negative and 
positive feedback loops are possible only in the phase domains where superposition holds and only when they are accessed by 
nonlinear phase transitions. 

The small-signal linear dynamic range revealed in the observations gives us the possibility to decompose cortical impulse 
responses by fitting to them the sums of linear basis functions, which have been found to be solutions of ODE [2| and have 


learn from the finding of the wave-pulse duality in laboratory observation of brain dynamics is that we can 
our discoveries if we accept both the continuum of the waves and the discreteness of the pulses and use 
each determines the other: pulse and wave modes coexist simultaneously in circular causality. On the other 
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FIG. 8: Aa.n vs t for different values of n and “low” fc = 3. 


been shown in the present paper to be solutions of the couple of damped/amplified oscillator equations. Remarkably, the 
damped/amplified oscillator equations are at the basis of the formulation of the dissipative quantum model of brain 
Ell: EO, EJ] • It seems thus that the very same duality particle/wave — discreteness/continuum — built in in QFT upon which 
rests the many-body dissipative model opens a perspective in dealing with pulse-wave duality combining the complexity of 
nonlinear interactions with linear, stationary, ergodic solutions at the observational level. We stress that in this paper, as well 
as in the dissipative quantum model of brain, neurons and glia cells are considered to be classical objects. The quantum degrees 
of freedom are associated with the electrical dipole vibrational modes of water molecules and other molecules of biological 
interest. The observational level represents the macroscopic manifestation of the underlying dynamics at the elementary 
constituent level [35| . 

In this paper, the connection of the damped/amplified oscillators with the ODE formalism is made explicit by considering 
specific applications. According to our discussion, the duration of a particular neuronal (“memory”) process can be related 
to the values of the damping coefficients that describe the relevant neurological process, coefficients that can be determined 
from a plot of the poles of the fitting functions in s-space. Observations show that there are two mechanisms for determining 
lifetime. In the chaotic gas-like phase the lifetime is determined by the decay rate of the exponential impulse responses, which 
in turn is specified by the negative feedback gain. In the ordered liquid-like phase, the gain is specified by the distribution of 
the characteristic frequencies around the mean value of the carrier frequency (^. EtI]. 

We have that couples of damped/amplified oscillator solutions, {n, kj} and {n, kj}, are linked in a dissipation/absorption 
fashion mode. The key of this link lies in the common value of the parameters defining the system {Lk^ ,0Jkj,n}- Furthermore, 
adding modes in pairs can be easily represented in s-space. With infinitely many k modes, one can always pair them up. Then, 
if one kj mode, for a given j, has a shorter lifetime than the other one in the couple, the linked system becomes unpaired 
(“widowed”). However, since kj and kj, for each j, belong to continuum sets, there is also the possibility of creating new links 
(new couples) by selecting a new companion for the survived one. Therefore, such a possibility of recovering destroyed links 
could be explored. 

Remarkable is also the fact that couples are among modes represented in left and right side of the root loci, which thus 
includes the transitions traversing the imaginary axis. Actually, since kj (and kj), for each j, is a variable in a continuum 
range, then, we may also have dynamical couples of modes with different i and j, for any i and j, with continuously changeable 
i and j in their domains of variability. These would possibly describe an actively correlated neural population, as it happens in 
amplitude modulated neuronal sets observed in brain functioning [3, @, EH- In these neuronal assemblies, due to the continuous 
interchange of partners in the couples, nothing is ’’wired” to something else, we have a neuronal population made of unwired, 
constantly changing neuronal connections. This is truly a dynamical assembly of correlated neurons in laminated neuropil, 
which is unified by an electrical interaction known to neurophysiologists as ephapsis , and which we conceive as forming a 

continuous field of exchange of information. These assemblies are described by coherent states in the dissipative quantum model 
of brain (cf. [1,[3:E3l for mathematical details). Our experimental evidence in the material domain for this “superconducting¬ 
like” layer is the intermittent synchronization of oscillations over the entire scalp EEG [43| . and it gives us the immediate unity 
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life-time 



FIG. 9: Aao.n vs t for different values of n and “high” k = 30. 


of rich detail in feelings and insights that constitute perception. 

We finally observe that onr data offer insight into the mind-body problem. The mathematics of dissipative systems requires 
a double, which is a mirror image defined by that portion of the universe that is engaged by the subject in action. Our data 
suggest that the synaptic structure of memories goes forward in time, while the remembrances in the double fed by memories 
recede into the past. Both structures evolve independently, yet they are ineluctably entangled at each moment of creation. 
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APPENDICES 

A. Empirical data, pulse density and wave density 

The empirical setting and the scenario in which the data analyzed in this paper were collected are described in this Appendix. 
We will comment on the linear, stationary, ergodic dynamics of the activity function observed in either of two modes, wave 
amplitude of dendritic current or pulse frequency, and will consider the transfer dynamics from point attractor to limit cycle 
attractor. 


1. Linear, stationary, ergodic dynamics 

Consider a population of excitatory neurons interconnected with a population of inhibitory neurons. They are distributed 
in a plane forming an area of cortex. Each neuron has a branching dendritic tree providing a surface area for receiving input 
through synapses from « 10^ other neurons. Each neuron has a single axon tree with multiple branches to « 10"* other neurons. 
For single neurons we define an activity function that is readily observed in either of two modes, as wave amplitude of dendritic 
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current or a pulse frequency, depending on the mode of observation. For single cells the gold standard is a patch clamp on 
the surface or a microelectrode inside the cell to estimate wave amplitude of dendritic current intensity from measurements of 
postsynaptic potentials (EPSP and IPSP) or of the frequencies of axonal pulses (units, action potentials). There is an extensive 
literature (see e.g. 0 ) showing that within a self-regulated narrow range of cortical function the conversions between wave 
amplitude and pulse frequency are proportional and additive, i.e., linear, for pulse to wave at synapses and wave to pulse at 
trigger zones. 

Maintenance of that range is by a population of excitatory neurons sustaining background excitatory bias that linearizes the 
dynamics, giving the scale-free, temporal power spectral density, PSDt, 1//", that characterizes the EEG spectra. 

For single cells the difference is trivial. The modes have in common the real-time transmembrane voltage at the trigger 
zone where the axon emerges from the cell body. For populations the difference in modes is profound, because of the problems 
of measurement. Impulse responses are observed and measured by fitting them by nonlinear regression with sums of basis 
functions @|. For a cortical population we define a wave density function and a pulse density function, which we sample over 
their time and space distributions in and around cortex. The necessity for doing so inheres in the fine structure of cortex, the 
neuropil. The density of cells is extreme (10® neurons/mm^, area in humans of 10® mm^ and perhaps 10 times that in glial 
cells), all of which contribute axons, dendrites and protoplasmic filaments between the cell bodies. With due attention paid 
to the architectures of cells and layers and the geometry of sources and sinks of dendritic ionic currents, we estimate the wave 
density from the amplitude of the local field potential (LFP) or electroencephalogram (EEG), which is the spatial sum of the 
contributions from all of the active neurons in the local area. 

It is not possible with prevailing techniques to record simultaneously sufficient pulses (« lO'*) to estimate the pulse density 
function. Under the assumption of ergodicity and stationarity, we record long trains of pulses of single neurons, from which 
we estimate the relative frequency of firing by calculating the probability of firing conditional on the amplitude of the wave 
density @]. For test inputs we deliver single electric shocks to the axonal pathways to or from cortex (see Ref. Ql for details). 
The impulse input delivers all frequencies. The impulse responses (average evoked potentials (AEP) in the wave mode and 
post stimulus time histograms (PSTH)) in the pulse mode) display the characteristic frequencies of the cortex averaged over 
the time span required to collect a sample. The two modes of activity are separated for observation by band pass filtering and 
time and space ensemble averaging, so that the generic equation for the impulse response in time and space is the sum of the 
responses in the two modes of observation: 


0{t,x,y,p,v) = Op{t,x,y,v) + Ov{t,x,y,p), (Al) 

where p and v stand for pulse density and wave amplitudes, respectively. 

By paired shock testing at differing intervals and intensities we test for domains in which additivity and proportionality 
hold. It turns out that Q] test inputs that give impulse responses within the self-regulated amplitude range of the EEG reveal 
a small-signal linear dynamic range. Superposition enables us to decompose the cortical impulse responses by fitting to them 
the sums of linear basis functions, which are the solutions to linear ordinary differential equations (ODE). This finding opens a 
major portal into the understanding of the complicated dynamics of cortex, mainly because it facilitates analysis and modeling 
of the multiple feedback loops that we encounter in even the simplest forms of cortex. 


2. Activity is simultaneously pulse density and wave density 

With these properties the general equation to express the activity density for wave or pulse can be separated into sequential 
operators in time F{t), space H{x,y), and amplification or attenuation Gp{v) at axon trigger zones or Gv{p) at synapses @|: 

Op{t,x,y,v) = F{t)Gp{v) H{x,y), at trigger zones, (A2) 

Ov{t,x,y,p) = F{t) Gv{p) F[{x,y), at synapses. (A3) 

For simplification a lumped model is developed based on a representative sample taken from a point in the population, and 
the spatial dimensions are subsumed by the selection and fixation of the stimulating and recording electrodes, so F[(x,y) = 1. 
The intrinsically nonlinear pulse-to-wave conversion at synapses and the wave-to-pulse conversion at axonal trigger zones are 
replaced by coefficients, kij, that express the linearized gains in the small signal range. We find that in all instances where 
it is possible to measure impulse responses in the two modes, the form in both modes contains a damped cosine in which the 
frequency w, phase of onset (f, and decay rate a are equal, and the amplitudes v or p differ only by an arbitrary constant for 
conversion of units of measurement which is set to one: 

0{t) = Ov exp{—at) cos{ujt + (f) + Op exp{—at) cos{ujt + (f) (A4) 

The range is extended into the nonlinear domain by piece-wise linearization and display of root loci in the complex plane 
of the Laplacian operator as functions of the stimulus parameters or of biological control parameters that are imposed by 
behavioral, physiological, surgical and pharmacological manipulations of the brain Q|. 

The wave-to-pulse nonlinear gain function was discovered to be asymmetric. Thereafter the values of the forward gains, kij, 
in the function Gp{v) were determined by the slope of the tangent to the nonlinear gain curve at the operating point Q]. 
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3. Attractor dynamics: from point attractor to limit cycle attractor 

A weak electric shock to the main input pathway to the prepyriform cortex from the olfactory bulb, which is called the 
lateral olfactory tract (LOT), excites a small subset of excitatory neurons. These few excite a larger number of inhibitory 
neurons, which in negative feedback inhibit a still larger number of excitatory neurons. They dis-excite inhibitory neurons, 
which dis-inhibit excitatory neurons, initiating an impulse response with the form of a damped cosine (Eq. (El). The impulse 
responses of the excitatory and inhibitory populations are in quadrature, having the same characteristic frequency but with the 
excitatory population leading on average 7r/2 radians (90°). The exponential decay to the baseline is modeled by postulating a 
point attractor under perturbation by noise. Convergence to the pre-stimulus baseline is a necessary condition for observation 
of both modes. The impulse response must terminate before the next impulse is given. The decay rate, a, must be negative, 
and the characteristic frequency, uj, must be given by the conjugate pair of poles closest to the imaginary axis of the complex 
plane. 

The frequency and decay rate of the complex conjugate pair of poles is primarily determined by the negative feedback gain, 
kn, between the excitatory population, KI^, and the inhibitory population, KU. It is strongly modulated by the positive 
feedback gains, fcee and ku, within respectively the A7e and Kh populations. 

Modeling the dynamics of the three loops can be postponed by noting the fact that the neurons in the populations get no 
information about what are the sources of their inputs. They are embedded in neuropil and integrate whatever synaptic input 
they receive from all sources. This fact can be modeled by postulating that when a neuron fires a pulse, the axon distributes 
the pulse to « 10^ other neurons, and some of those to « 10® others, and so on, reverberating through the neuropil as it 
decays. The neural activity that feeds back to the originating neuron (or neurons) embedded in the neuropil can be modeled 
as a one-dimensional diffusion process, using the linear transcendental Laplacian operator, K 2 exp[—{sT 2 )^'^], where K 2 is the 
lumped gain and T 2 is the lumped time delay of the distributed feedback. 

The advantage of this simplified model was that it revealed a fundamental property of the interactive populations of excitatory 
and inhibitory neurons, modeled as the Kllsi set. When the impulse input intensity was fixed in resting subjects, repeated 
samples of 0{t) revealed spontaneous variation of the damped cosine. It turned out that K 2 and T 2 were linearly correlated. 
Variation in lumped gain K 2 gave root loci for lj and a that crossed into the right half of the complex plane, reversing the sign 
of a from attenuating to amplifying, and predicting transfer of cortical control from a point attractor to a limit cycle attractor. 
In physiological terms the closer cortical axons were brought by noise to their thresholds, the more likely they were to fire. 
This well-known biophysical nonlinearity at the microscopic level appeared at the macroscopic level as the root locus in Mode 
2 [2], thereby revealing a mechanism capable of explaining how a surge in neural activity could provide the transition energy 
required to initiate sustained oscillation by a phase transition resembling a Hopf bifurcation. 


B. Explicit examples of computations 

We present explicit computations in terms of Eqs. and m for two examples from Ref. Q], Section 2.5.3: the average evoked 
potential from KII set and poststimulus time histograms of single glomerular neurons in olfactory bulb. We also summarize 
the discussion on the existence of the limit cycle. Preliminary results have been presented in |45l| . 


1. Average evoked potential from KII set 


We refer to Eq. ([1611 and will drop the kj suffix for convenience, keeping only the label of each momentum in order to repeat 
the analysis as in Q], assuming that the /3’s are common for each mode (i.e. j3i = ft). 

The prepyriform response to a 5{t) input can be modeled by 2 parallel filters where each is cascaded into a lead-lag filter [^. 
This means that the response function c{t) = ci{t) -|- C 2 {t) has one real zero, one real pole and two complex conjugates poles 
per component, which can be written in the s space as 


C’(s) = Ci(s)-bC'2(s) , 


Kj{s + Zj) 


(s + B,)(s + B*)(s + ft’ 


j = 1, 2 , 


(Bl) 


where the K’s are constants, the 7’s can be found by graphical analysis and determine the phase of the output. By partial 
fraction expansion, Eq. dnn), when separated in modulus and phase, yields 




+ ■ 


Pi L rf + (/? - 71)2 
ri(^i-ft 


]’( 


2(s + Bl) 2(s + B^) 

)■ 


-b 


[Pf + (^1 - 7i) 2] i [Pf + (/? - 7 i) 2] i (s + ft 

with 0i = tan~^{Im{Ci{s))/Re{Ci{s))). The same result applies to <72(5). 


(B2) 
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Since we are dealing with a S{t) input and c{t) output, we have the input V'(s) 
gives 


C{s) so that use of the Laplace transform 


Ci{s) 


l/i(s) = Vi[ 


s + Bi 


+ 


s + 



(B3) 


and similarly for C 2 {s). This allows us to compare the experimental data described by Eq. (I16II with the modeled output 
(Eq. (IBll) i and find the values of Comparison of Eq. (IB2I) and Eq. (IB3I) shows that: 


Ki = VlEi 

sin(vJi) = ■ 


rf + (/3 - 7i)^ 1 ^ 

rf + (^i-7i)2j ’ 


01 = -(^1 - -) 


[r? + (2i-7i)^]^[r? + (/?-7in^ ■ 

with an equivalent expression for the if2-term. Then, the full expression of C(s) can be written as 


, _ {Ki + K 2 )s^ + Pis^ + P2S + B3 

(s + B\)[s + BJ)(s + B2){s + )('S T P) 

with Pi = A'i(zi + 272) + 1^2(22 + 271) , 

P2 = P'1 (27221 + 72 + ri) + P2(27122 + 7i + rf), 

P3 = Pl2l(7|+r2)+P222(7?+r?) . 


(B4) 

(B5) 


(B6) 

(B7) 

(B8) 

(B9) 


Zeros and poles of Eq. (lB6t in the complex plane are plotted in Fig. 2. We can regard each damping/sinusoidal-damping term 
in Eq. (Hi) as the observable excitations of a particular fc-mode, being the modeled output a linear superposition of different 
fc-modes. Moreover, we will postulate the double of this system as an amplihed version of Eq. Hi : 


v{t) = Vki [sin(rfcit + - sin(93fci)e + 

+ Vfc2[sin(rfe2t + - sin(v2fe2)e“^''2‘] , (BIO) 


where the Pkj are the same for the amplihed and the damped system. 

Each coefficient will satisfy the restrictions given by Eqs. mi) - (HSJ. Moreover, by requiring that each /3 value is a solution of 
the same fc-mode as the sinusoidal damping/amplifying term, when substituting the Pkj term into Eqs. (H)) and (HI we obtain: 


Pkj ,n — 7fej ^TT, , 


Pkj^n ^kj,n^^^kj,n * 


(Bll) 


Thus, in order for j3kj to be common to (i.e. the link between) the damped/amplihed pair, it has to be 'yki,n = —7^^ „ and 
which agrees with the previous discussion. Therefore, /3 (which is always a damping coefficient) is indeed a link 
between the damped/amplihed coupled mode. 

If one substitutes each value of /3 for a single kj mode in Eq. (II6II , (dropping the k and n subscripts again for simplicity) the 
solutions reduce to a single damped/amplihed sinusoidal term with a phase factor: 


V{t) = sin(rt)e"'"‘e=^'"’ , V{t) = sin(rt)e'"‘e±'"’ . (B 12 ) 

We can now repeat the analysis of Ref. Q] for the addition of two different fe-modes but with an output of the form of 
Eqs. (IBI2II . Note that, according to the discussion in @], the absence of a decaying j 3 coefficient implies the cancellation of 
the real zero and the real pole in the transmission function of the output in s-space (i.e. z = P for each fc-mode) so that the 
model simply reduces to two parallel band-pass hlters which are not cascaded into a lead-lag filter, represented by 2 complex 
conjugates poles in s-space. The modeled output c{t) = Ckj^{t) + Ck^{t) leads hnally to 

C'(s) = C'i(s)-fC2(s) = 

^ (lEi -b K 2 )s^ + (R'i 272 -b Jf227i)s -b + Tj) -f Jf2(7? + Ef) 

(s -I- 71 -I- iEi)(s -I- 71 - zEi)(s -I- 72 + *E2 )(s -I- 72 - fE2) 

in s-space. The numerator of Eq. (IB 13 | is a quadratic equation with two complex conjugate roots 
resemblance with Eqs. m and by using the dehnitions of L and lo: 

{Ki + K 2 )s^ + (iEiL2 + K 2 Li)s + {Kioji + K 2 ujI) = 0 . 

Therefore, we conclude that a linear superposition of two fc-modes of the form of Eqs. (HSI, where each satisfies Eq. © with 
a different value of /?, leads to 4 complex poles and 2 complex zeros in s-space. Then, each extra mode of the same type added 
to the output function will increase the number of complex poles and complex zeros by two when the output is represented in 
s-space. In addition, when the amplified terms are also considered (solutions of Eq. ®), its poles and zeros are a reflection in 
the imaginary axis of the ones presented by Eq. (IB 13 I) . 


(B13) 

. In fact, one can notice the 
(B14) 
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2. Poststimulus time histograms of single glomerular neurons in olfactory bulb 


The second example we consider is the modeling of the synapses of axons in the primary olfactory nerve (PON) (cf. also 

[3], Section 5.2.2). 

Starting from a model based on a characteristic set of differential equations, the transfer function for the feedback subset of 
the closed loop representing these synapses is 


_ KcKe8.5 X 10^*^(2300 -3)^_ 

(s + 7 + ir)(s + 7 — fr)(s +pi)(s +P 2 )(s + P3)(s +P 4 ) 


(B15) 


where all the coefficients are real and positive, Kc is the closed loop gain, Ke is the feedback gain and the pole coefficients are 
experimentally determined. The inverse transform of Eq. (IB15II is 

Pc(t) = + Pse”^^* + + Ps sin(rt + (p)e-'"‘] , (B16) 

where the P’s are parameters to be determined and Kg is the overall gain. The validity of Eq. (IB15II is confirmed by the 
similarity of Eq. (IB16I) with the recorded poststimulus time (PST) histograms. The overall transfer function is 


P,{s) 


_ KgK^8.5 X 10^'^(s + gi)(s + 22)(2300 -s)^_ 

(s + 7 + ir)(s + 7 — ir)(s +Pl)(s +P 2 )(s +P3)(s +P4)(s +P5)(s +P 6 ) 


(B17) 


where zi,Z 2 = -(P 5 +^ 6 )±d ^ j-ggpgg^^jygjy^ with d = [(ps + pe)^ — 4p6P6(l + K'}] ^ ■ K' is the amplitude of the surge relative to 
the pulse. Note that we can calculate the zeros of Eq. (IB17I) from 21,22 and the poles will change with so that we can 
compute the root locus plot. Again, the transfer function is 

Pq{t) = Kg[Pie~^^^ + P 2 e“*"^* + Pze~^^^ + P 4 e“^"‘‘ + P 5 sin(rt + <p)e~'^^+ 

+ Pee"^'^ + Pye"^*'] , (B18) 


which is very similar to the recorded PST histograms. 

The transfer functions in t-space are then interpreted as a linear superposition of the nonlinear fc-modes. The position of 
the poles of the transfer function in s-space which are gain-dependent determine the coefficients in the superposition (Fig. 4). 


3. Limit cycle discussion 


The transfer function D{s) for a Kill set output from a KII set with negative feedback from a KI set is given by: 


D{s) = 


C{s) 

lpK,C{s)P{s) ’ 


(B19) 


where C(s) = (6.25 x 10®)/(s-I-«250)(s — i250)(s-f 100) and P(s) = 100/s(s-I-100) are, respectively, the KII and KI set 
transfer functions. As Ke changes, the poles of Eq. ()B19I) change accordingly: 


D(s) 


_6.25 X 10'^s(s -f 100)_ 

(S-|-pi)(s-|-P2)(S -|-P3)(S-|-P4)(S-bps) 


(B20) 


where the p„ poles can be real or complex. The D{s) root locus is plotted in Fig. 5. The larger poles cross the imaginary 
axis from the left to the right as decreases. This signal an unstable limit cycle. At frequencies corresponding to purely 
imaginary poles (« 160 rad/sec) the limit cycle is stable. 
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